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Abstract In this paper, we use duality arguments "a la Michel Pierre" to establish global 
existence of classic solutions for a class of parabolic reaction-diffusion systems modeling, for 
instance, the evolution of reversible chemical reactions. 

1 Introduction 



This paper is motivated by the general question of global existence in time of solutions to 
the following reaction-diffusion system 



(5) 



ut — diAu 
vt - d2Av 
wt — d^Aw 

du dv 
-{t,x) = -it,x) 

u{0, x) 
f (0, x) 
w{0, x) 



where is a bounded regular open subset of 



(0,+cxd) X n, 

(0, +cxd) X n, 



dw , , 



-w^ + vP'vf^ (0, +C50) X 0, 
(0,+oo)x5O, 



(^l) 

(^2) 
(^3) 



vq{x) > 
wq^x) > 



{di,d2,d3,a,l3,j) G (0,+(X))3 x [l,+oo)3. 



Note that the system (S) satisfies two main properties, namely : 
(-P) the nonnegativity of solutions of (5) is preserved for all time ; 

(M) the total mass of the components u, v, w is a priori bounded on all finite intervals (0, i) 



If a, f3 and 7 are positive integers, system (5) is intended to describe for example the evolution 
of a reversible chemical reaction of type 

aU + pV ^ -fW 

where u, v, w stand for the density of U, V and W respectively. 

This chemical reaction is typical of general reversible reactions and contains the major difficulties 
encountered in a large class of similar problems as regards global existence of solutions. 
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Let us make precise what we mean by solution. 

By classical solution to (5) on Qt = (0, T) x il, we mean that, at least 

(i) {u,v,w) eC{[o,Ty,L^nf)nL°°{[o,T] X nf,yT e (o,r); 

(ii) yk,£= l...iV, VpG (l,+oo) 

dtu, dtv, dtw, dx^u, d^^v, d^^w, dx^x^u, dx^x^v, dx,^xeW, u,v,w € LP{{0, T) x n) ; 
(Hi) equations in (S) are satisfied a.e (almost everywhere). 

By weak solution to (S) on Qj- = (0, T) x J7, we essentially mean solution in the sense of 
distributions or, equivalently here, solution in the sens of the variation of constants formula with 
the corresponding semigroups. More precisely 



u{t) = SdAt)uo+ SdAt-s){w''{s)-u''{sy{s))ds 
Jo 

vit) = Sd,{t)vo+ [ Sd,{t-s)iwUs)-u"{sy{s))ds 
Jo 

w{t) = Sd,{t)uo+ [ Sd,{t-s){-w^s) + u''{s)v^{s))ds 



where Sd^^-) is the semigroup generated in L^(Q) by — djA with homogeneous Neumann boun- 
dary condition, 1 < i < 3. 

By just integrating the sum (£'i) + (£'2) + ^{Es) in space and time, and taking into account 
the boundary conditions A{diu + d2V + d^w) = 0^ , we obtain 



u{t) + v{t) + 2w{t) = uq + vo + 2wo t>0. (1) 
1 Jn 

Together with the nonnegativity oi u, v and w, estimate ([1]) implies that 

Vt > , \\u{t)\\Li(n), \\v{t)\\L^(n), lk(0llLi(f7) < ||no + wo + 2u;o||Li(n)- (2) 

In other words, the total mass of three components does not blow up; u{t), v{t) and w{t) rest 
bounded in L^[Q) uniformly in time. 

Although one has uniform L^-bound in time, classical solutions may not globally exist for 
diffusion coefficients di, d2, d^ which are not equal (global existence obviously holds if di = d2 = 
ds). As surprisingly proved in [12j and [16], it may indeed happen that, under assumptions (P) 
and (M), solutions blow up in finite time in ! In particular, classical bounded solutions do 
not exist globally in time. 

If uo,vo,wo G L°°(ri), local existence and uniqueness of nonnegative and uniformly bounded 
solution to (5) are known (see e.g. [I3]). More precisely, there exists T > and a unique classical 
solution {u,v,w) of (5) on [0, T). If Tmax denotes the greatest of these T's, then 

(Tmax < +00) =^ lim (ll^i(t)llioo(Q) + \\v{t)\\Loo(_Q) + \\w{t)\\ loo = +00. (3) 

t/ J- max 

To prove global existence {i.e. Tmax = +00), it is sufficient to obtain an a priori estimate of the 
form 

VtG[0,rma.), \\u{t)\\L^^a) + \\v{t)\\L^^tt) + \\w{t)\\L^(^^)<H{t), (4) 
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where H : [0, +cx)) [0, +00) is a nondecreasing and continuous function. 



This type of estimates is far of being obvious for our system except the case where diffusion 

u + V + 2w satisfies 



coefficients di, 


^2, 


are equal i.e di = d2 


= d3 = 


d. Indeed, Z = 








r Zt - dAZ 


= 


(0,+oo) X 






{E)< 


dZ 

dn 
[ Z{0,x) 


= 
= Zo 


(0,+oo) X 
[x) X £ il, 


where Zq{x) = 


uq{x) 


+ vo{x) 


+ 2wo{x). 







In particular, we deduce by maximum principle that 

\\u{t) + v{t) + 2?i;(t)||i«=(f^) < ||no + ■uo + 2t(;o||L°o(n), t > 0. 
Together with nonnegativity, this implies 

\\u{i)\\L^{Q) + \\v{t)\\L^{Q) + \\w{t)\\L^(^Q) < \\uo + vo + 2u;o||l=c(q), t > 0. 

In other words, u{t), v{t) and w{t) stay uniformly bounded in L°°(il) and therefore T^ax = +00. 

In the case where the diffusion coefficients are different from each other, global existence is 
considerably more complicated. It has been studied by several authors in the following cases. 
First case a = /3 = 7 = 1. 

In this case, global existence of classical solutions has been obtained by Rothe [13] for di- 
mension N < 5. Later, it has first been proved by Pierre |lUj for all dimensions N and then by 
Morgan [9]. 

The exponentiel decay towards equilibrium has been studied by Desvillettes-Fellner in 
the case of one space dimension. 

The global existence of weak solutions has been proved by Laamri [7j for initial data uq, vq 
and wq only in L^($7). 
Second case 7=1 regardless of a and /3. 

In this case, global existence of classical solutions has been obtained by Feng [4j in all 
dimensions and more general boundary conditions. 
Third case a + /3<2or7<2. 

In this case, Pierre [11] has proved global existence of weak solutions for initial data uq, vq 
and Wq only in L'^(Q). 

Our paper mainly completes the investigations of [|1], |9], [10], [13]] and [[7], [H]]. As far as 

we know, our results are new either when a + /3 < 7, or when 1 < 7 < — — — regardless of a and 

' iV + 2 

(3. For the sake of clarity, we decided to focus in this work on the question of global existence in 
time of solutions in the case of homogeneous Neumann boundary conditions. So, we shall prove 
global existence of classical solutions to system (S) in the following cases : 

* a + /3 < 7; 

* {di = d^ or ^2 = da) and for any (a, (3, 7) ; 

* di = d2 and for any (a, /3, 7) such that a + /3 7^ 7 ; 

1 < 7 < ^2 

For the sake of completeness and for the reader's convenience, we shall also give a direct 
proof different from that of Feng [2] in the special case 7 = 1. 
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Notation : Throughout this study, we denote by Q's various positive numbers depending only 
on the data and for p € [1, +00 [ 



i/p 



u{t,x)\Pdx\ , |in||ip(Q^) 



u{t, x)|^ dtdx I , 



T 

/ \u{t,x)\^ dtdx 
n J \Jo Jn J 

u{t)\\oo = esse sup2.eo|'u(t,x)[, \\u\\ ^oo (^q^^ = esse sup(j^^.)gQ^|n(t,2;)|. 



2 The main results 

One of the main ingredients for the proof of our results is the following lemma which is 
based on the regularizing effects of the heat equation. This lemma has been introduced by 
Hollis-Martin-Pierre in [5]. 

Lemma 1 Let T > and {(f), ip) the classical solution of 



(j)t - diA(l) 
tpt - d2A'ip 

cPiO,x) 
[ V'(0,x) 



{t,x)e{o,T)xn 
g{<p,i>) (t,x) e (o,r) X o 

G (o,r) X dVt 






i}Q{x) 



{t,x) G (o,r) X dn 

X €Q 
X €Q. 



Assume that f + g = 0, then for each p € (1, +00), there exists C such that for all t G (0,T) 

MiPiQ,) < c [UhnQt) + 1] ■ (5) 

A more general version of this lemma can be founded in [IH lemma 3.4]. □ 
2.1 The case a + /3 < 7 

Theorem 1 Assume that < uo,vo,wo < M where M is a positive real. 
If a + f3 < ^, then the system (5) admits a global classical solution. 

Proof : 

• Let T G (0, Tmax) and let t G (0,T]. Thanks to the nonnegativity of u, v and w, we deduce 
from the equation (Ei) that u is bounded from above by the solution U of 



(^1) 



f Ut-diAU 
I UiO,x) 



w"/ {t,x) e {0,T) X n 

{t,x) e {0,T) X dn 

uo{x) X G Cl, 

and we deduce from the equation {E2) that v is bounded from above by the solution V of 

w"/ it,x) e iO,T) xQ 

G (o,r) X 

Vo{x) X G 0. 



r vt-d2Av 

dV , , 

on 
I ^(0,x) 
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Therefore it is sufficient to show that w E U'{Qt) for p large enough. 

• Let g > 1. Multiplying the equation (£^3) by w'' and integrating over Qt, we get 

^— / w'^+\T)+qd3 [ f \Vw\^w''-^+ [ f w^^^ = f [ u'^v^w'i + Ko (6) 
9 + 1 J Jqt J Jqt J JQt 

where 



Thanks to Holder's inequality, we have 



Qt 
where 

- + - + — ^ = 1. 

r s 9 + 7 

Since a + /3 < 7, we can choose r such that ra < q + j and s such that s/3 < <? + 7. To convince 

7 

oneself, it is enough to draw the straight line with cartesian equation x + y = and to 

Q + 1 

a (3 
identify the points with coordinates ( ,0) and (0, ). 

q+1 q+i 

Then L5+'>'(Qt) C L'^'^Qt) and L1+^{Qt) C L^'iQr)- Consequently, there exists Ci such that 

uVw;" < Ci||^x||2,+,(Q^)|b||J,+,(Q^)||u;||^,+,(Q^). (8) 

By virtue of lemma 1, there exists C2 such that 

Iklli^+nQT) ^ '^2(1 + \\w\\li+i{Qt)) (9) 

and there exists C3 such that 

\Ml-i+i{Qt) ^ '^3(1 + \\w\\l-i+i{Qt))- (10) 
Thanks to Q and (1101) . estimate dH]) can be written 



V^V^W"^ < C4 (1 + ||w||i7+9(Qy))" (1 + \\w\\l~<+i{Qt)) (1 + \\w\\l~<+i{Qt))'^ ■ (11) 

If ||il'||l7+9(Qj,) < 1 then the proof ends up. Otherwise, there exists C5 such that 

n-v^w''<C,\\w\\lt,:t^^^y (12) 



So we deduce from 



w'^+^<C,\\w\\l+::i^^^+Ko. (13) 



Q 



T 



With the notation R := j I w'^'^'^ , estimate can be written 

'Qt 



R<C5R^^+Ko. (14) 
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Since q + a + l3<q + ^, by applying Young's inequality to p^ . we obtain 

(l-e)i?<i^o + C6. (15) 
Then, for e G (0, 1), we have the desired estimate 

\\M\L'i+i(QT) ^ C'7- (16) 

Going back to (Pi) and (-P2)) we have, by choosing q such that ^ > — l!l — and thanks to 
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the L^-regularity theory for the heat operator (see [8j), 

||m||l°°(Qt) ^ C's (17) 

|b||L«=(QT) < Cg. (18) 

Now going back to (-E3), we deduce from (fT7|) and (fT8|) that there exists Cio such that 

||w||ioo(Qj,) < Cio- (19) 

This implies that Tmax = +00. □ 

Remark This method seems to be specific to the case a + /3 < 7. It fails when a + /3 > 7 since 
some restrictions on the parameters q, /3, 7 and on the diffusion coefficients will appear. 

2.2 Case where di = or d2 = d^ or di = ^2- 

Theorem 2 Assume that < uq^vq,wq < M. 

(i) If di = d-i or d2 = d^, then system (S) admits a global classical solution for any (a, f3,j). 
(a) If di = d2, then the system (S) admits a global classical solution for any (a,/3,7) such that 

a + /3/7- 
Proof : 

(i) Assume that di = ^3 = d, we have 

{u + w)t - dA{u + w) = 0; = ; (u + w){0, x) = no(x) + u'o(x). 

on 

We deduce by maximum principle 

\\u{t) +w{t)\\ao <\\uo + Wo\\oo- (20) 

Together with the nonnegativity of u et lu, this implies that u{t) and w{t) are uniformly bounded 
in L~(f]). 

By going back to (£'2) and thanks to the L^-regularity theory for the heat operator (see [6]), we 
conclude that 11 

^(^)||cxD IS unifornily boundGd. in -L^(r2) on a,ll interval [O,?^] so tlia-t 2^max 

= +00. 

(a) Assume that di = d2 = d. The case a + /3 < 7 was already handled in the theorem 1, so 
it remains only to tackle the case 7 < a + f3. Moreover, one can assume that uq 7^ vq since if 
uq = Vq the result is obvious. 
Since di = d2 = d, we have 

diu — f) 

(n — v)t — d/S.iu — w) = ; = ; (n — t')(0, x) = u^ix) — voix). 

on 
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The maximum principle then imphes \\u{t) — v{t)\\oo < \\uq — foUoo = C. Hence we have 

= u'^v'^ + n"/3(en + (1 - e)vf~\u - v) where 9 g]0, 1[ 

Thanks to Young's inequahty, there exists Cn > and C12 > such that 

C7iiu"+^ < n^i;^ + C12. (21) 
By virtue of (PT|) . equation (E'l) imphes that 

ut - diAu + Cuu"^'^ < w'' + C12. (22) 
Let q > 1. Multiplying ([22]) by u'^ and integrating over Qr, we obtain 

^ /" t.^+i(7^)+'7c^2 / / [V^xpn^-i+Cii /" / n^+"+^ < [ [ w^u'^+Cu [ [ W^+Ki 

"T i Jfi J JQt J JQt J JQt J JQt 

(23) 

,9+1 



q + 
where 



9 + 

Thanks to Holder's inequality, we have 



tw^n" < ( / / i W^') (24) 



a + /3 + g , Q; + /3 + g 

where r = and s 



7 q + a + f5 - J 

Lemma 1 implies that there exists C13 such that 



l/r 



^7 



If ||ti||L<7+a+/3(-Qj,) < 1 then the proof ends up. Otherwise, there exists C14 such that 

l/r 



w^^j <C^,Ml,^.^,s^Q^y (25) 
Since qs < q + a + (3, we have L'^"'""^'^((5t) C L'^^{Qt), then there exists C15 such that 

N I/S 

^^-^^J <CMU^^,^Q^y (26) 
Denote 5 := / / n«+"+^. Estimates ([25]) and ([26]) imply that 

w^u'i <CieS^^. (27) 

Qt 



Moreover, since 

L1+'^+P{Qt) c L1{Qt), there exists Cn such that 

Ci2 / / -u*? < CiT^^™. (28) 



Since 7 < a + /3, by applying Young's inequahty, there exists Cis such that 

Cie^™ < |5 + Ci8. (29) 

Applying again Young's inequality, there exists C19 such that 

C7i75^+^ < |5 + Ci9. (30) 

Consequently, estimate ([23]) implies 

{Cii-e)S<Cis + Ci^ + Ki. (31) 
By choosing e < Cn in (f3T|) . there exists C20 such that 

ll'"llL<?+"+/3(Qr) < C2Q. (32) 

Thanks to lemma 1 and estimate (|32p there exists C21 such that 

ll'"^llL9+«+/3(Qr) < C21- (33) 

By going back to (Pi) and {P2), we have by choosing q such that — ^ > — — and thanks 
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to the L'^-regularity theory for the heat operator (see [6]) 

\\u\\l^(Qt,) < C22 (34) 

\\v\\l^(Qj.) < C23- (35) 

Now let's go back to {E3), we deduce from and ([35]) that there exists C24 such that 

lkllL°°(QT) ^ ^24- (36) 

This implies that T^^x = +00. □ 

Remark Even in the last case i.e di = d2, global existence or blow-up in the limit case a-|-/3 = 7 
remain an open problem. □ 

2.3 Case 1 < 7 < regardless of a and B. 

N + Q 

Theorem 3 Assume that < uq,vo,Wo < M where Af > 0. Ifl<j< ^ ^ ^ , then the system 
(S) admits a global classical solution for any {a, 13) G [l,-|-oo)^. 
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Proof : 

Let T € (0, Tjnax) and let t € (0, T]. Thanks to the nonnegativity of u, v and w, we deduce from 
the equation (Ei) that u is bounded from above by the solution U of 



f Ut-diAU 
dU , , 



w 



{t,x) G (o,r) X n 

{t,x) e iO,T) X dn 
uo{x) X G Q. 



Therefore it is sufficient to show that w € L^{Qt) for p large enough. For this we have to 
distinguish the case 7 = 1 and the case 7 > 1. 
• Case 7 = 1 and a, /3 > 1. 

Let us recall that global existence of classical solutions for {S) when a = /3 = 7 = 1 has been 
studied by several authors. It has been obtained by Rothe [13j for dimension N < 5. Later, it 
has first been proved by Pierre [10] for all dimensions and then by Morgan [9j. 
Independantly, Feng [4J has proved global existence in the case 7 = 1 regardless of a and (3 and 
more general boundary conditions. 

For the sake of completeness and for the reader's convenience, we give here a simple and direct 
proof in the last case (7 = 1 regardless of a and In our proof, we use an idea introduced by 
Pierre in PD] and applied in [S]. 



For any p > 1, we deduce from (Pi) and the semigroup property 



\u{t)\\p < \\uq\\p + \\w{s)\\pds. 



(37) 



By applying Holder's inequality for p > 1 and thanks to ([5]), we obtain 

i/p 

\\w{s)\\pds <t^/p' ( / / wPdsdx] 



,ds < t^'P' (^j^ j ujPdsdx^ < t^/P'C25 



1 + 



dsdx 



Jn 



i/p 



(38) 



where p' 



P 



p — 1 



h{s) ds 



For t € (0, r], let us set h{t) := / \u{t,x)\P dx. Inequality ([57]) can be written 

Jn 

Taking the p^^ power of (|39p we obtain 

h{t) < 2P'^CIq + 2P-^Clj f h{s) ds. 

Jo 

But, inequality ([I0|) is a linear Gronwall's inequality, then 

||u||iP(Qy) < C28. 

Repeating the method above with v instead of u, we obtain 

\\v\\lp{Qt) ^ ^29- 



(39) 

(40) 

(41) 
(42) 
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Estimates (|41|) and (j42|) imply that for some q > — ^- — 

||n"7;^|U,(Q^) <C3o. (43) 
Going back to equation (£^3) we have, thanks to the L'^-regularity theory for the heat operator, 

||w^||l-(Qt) ^ C'si- (44) 
This concludes the proof for the case 7 = 1 regardless of a and /3. □ 
.Casel<7<^. 

The proof in this case is based on lemma 1 and these two following lemmas. 
Lemma 2 (Michel Pierre) Let T > and let Z the solution of 

( Zt-A{A{t,x)Z) < (t,x) G (0,r) X 0, 

dZ 

— (i,x) = {t,x) e {o,T) X dn, 

on 

[ Z{0,x) = Zo{x) xen. 

Assume that < d < A{t, x) < D where {d, D) € (0, +00)^. Then, there exists C = C{T, d, D, Q) 
such that 

||^||l2(Qt) - C\\^o\\L2{n)- 
For a general version of this lemma, see |1H proposition 6.1] or [3l theorem 3.1]. □ 

Lemma 3 Let {p,q) such that 1 < p < q < +00, d > and Sd{t) the semigroup generated in 
LP{Vt) by —dA with homogeneous Neumann boundary condition. Then 

-JV/l In 

\\Sdit)Y\\g < {C{n)m{t))~^p V \\Y\\p, for all Y £ LP{n), t > (45) 
where m{t) = min(l,t). 

For a proof of this lemma see for instance \i?>\ Lemma 3, p. 25] or Theorem 3.2.9, p. 90]. □ 
We now go back to the proof of theorem 3. 

By applying lemma 2 to the system (5) where Z = u + v + 2w and A = — ^^3^^ 

u + V + 2w 

we have u, v, w G L'^{Qt)- More precisely, there exists C32 such that 

IWWmQT)^ IHl^Qt)^ IMl^Qt) ^ C'32- (46) 

Now, we have thanks to the estimate (|45p with p > 1 and q = +00 

||n(t)||oo < \\uo\\o^ + C33 {t-s)^p\\w^s)\\pds. (47) 







By applying Holder's inequality, we obtain 



J {t-s)^\\w{sy\\pds< U {t-s)~ dsj U Ww^fiaWpds] . (48) 
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/■* -N N + 2 

We first remark that the integral {t — s)^^p~'^^ds converges when p > and we have 

Jo 2 

Jo Jo 







On the other hand, lemma 1 implies that 



(^1^ W{s)\\Pds^ = < C], (1 + ||u|Up.(q,))^ ■ (49) 

If ||'u||/^P7(Qt) < 1 then the proof ends up. Otherwise there exists C35 such that 

i/p 



WwHsW^ds^ <C-Ml,,^Q^y (50) 



Since 



/ r r \1/P 



'T 



it follows that (|47p can be written 



i/p 

uP^-P+' ] . (51) 



T 



If p(7 — 1) < 2, by choosing e € (0, min(p, l — pi^j — 1)), we deduce from (jl6|) and ([5T]) that there 
exists C37 such that 

||u||ioo(Qy) < C37. (52) 
2 4 iV + 6 

Note that the above condition 0(7 — 1) < 2 holds if7<lH — <1 + — = — . 

We establish in the same way that there exists C38 such that 

^ ^38- (53) 
Finally, for (£'3), we deduce from (j52p and (j53p that there exists C39 such that 

\\w\\l°-{Qt) ^ C'39- (54) 
N + 6 

This concludes the proof in the case 1 < 7 < . □ 

^ ' N + 2 

Remark : Our conjecture is that 7* = — — - is not optimal. In fact, when = 1 one can 
prove that the result of theorem 3 still holds for 7* = 7/2. □ 
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3 Conclusion 



• All our results are still true if we replace homogeneous Neumann boundary conditions by 
homogeneous Dirichlet boundary conditions, it suffices to replace lemma 3 by the following one. 



Lemma 4 Let {p,q) such that 1 < p < q < +oo, d > and S(i{t) the semigroup generated in 
LP{VL) by —dA with homogeneous Dirichlet boundary. Then 

-N, 1 1 ^ 

\\Sd{t)Y\\q < {ATTt)~^p^^' \\Y\\p, for all Y G LP{n), t > 0. (55) 

For a proof of this lemma, see for instance |14^ Proposition 48.4, p. 441]. □ 

• In the case where the diffusion coefficients are not equal (i.e. di ^ dj for all 1 < i 7^ j < 3), 
global existence of classical solutions for (5) or blow-up is still an open question when 

+ 6 

Our guess is that system (S) admits a classical global solution for all ^ ^ ^ < 7 < a + /3 and 
that there is a finite time blow-up when 7 = q + /3 and the dimension is large. ■ 
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